Phases of massive scalar field collapse 
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We study critical behavior in the collapse of massive spherically symmetric scalar fields. We ob- 
serve two distinct types of phase transition at the threshold of black hole formation. Type II phase 
transitions occur when the radial extent (A) of the initial pulse is less than the Compton wavelength 
{fJ.~^) of the scalar field. The critical solution is that found by Choptuik in the collapse of massless 
scalar fields. Type I phase transitions, where the black hole formation turns on at finite mass, occur 
when Xfj, ^ 1. The critical solutions are unstable soliton stars with masses < 0.6/i~^. Our results in 
combination with those obtained for the collapse of a Yang-Mills field [M. W. Choptuik, T. Chmaj, 
and P. Bizon, Phys. Rev. Lett. 77, 424 (1996)] suggest that unstable, confined solutions to the 
Einstein-matter equations may be relevant to the critical point of other matter models. 



PACS number(s): 04.25.Dm, 04.40.-b, 04.70.Bw 

The discovery of critical point behavior in gravitational 
collapse has highlighted the role played by non-linear dy- 
namics at the threshold of black hole formation, and has 
opened up a fascinating area of research in General Rel- 
ativity. 

Choptuik [|l[ performed the first definitive numerical 
study of critical behavior in the collapse of spherically- 
symmetric distributions of massless scalar field. His re- 
sults indicated that one parameter families of interpolat- 
ing solutions S[p] generically have a critical value p — p* 
such that (i) S[p < p*] are solutions in which the scalar 
field disperses to infinity, and (ii) S[p > p*] are solu- 
tions in which the field collapses to form a black hole. 
In slightly super-critical evolutions, Choptuik found that 
the black-hole mass has a simple power-law form 



Mbh -K\p-p* 



(1) 



where the critical exponent is 7 ~ 0.37 and if is a family 
dependent constant. For near critical evolutions the field 
asymptotically approaches a discretely self-similar form, 
with an echoing period A ~ 3.44 which is the same for 
all families, before either dispersing to infinity or form- 
ing a black hole. Based on these observations, Choptuik 
conjectured that a unique solution to the Einstein-scalar 
field equations acts as the intermediate attractor for near 
critical evolutions. Gundlach has directly constructed 
this critical solution, and has computed the echoing pe- 
riod and the critical exponent to be A = 3.4453 ± 0.0005 
and 7 = 0.374 ± 0.001 ||], thus confir ming the numer- 
ical estimates. It has also been argued that this picture 
is stable against the introduction of a small scalar field 
mass ^-|^. 

Motivated by the results of Choptuik, critical point 
behavior has also been studied in other models of gravi- 
tational collapse 1^1^. Dynamical self-similarity (either 
discrete or continuous) in near critical evolutions, and a 
scaling relation for black-hole mass, as in Eq. (pi), are 



common features of these models, although the numeri- 
cal value of the critical exponent 7 is model dependent. 
It is now well established that the power law form for the 
black-hole mass derives from the existence of a single un- 
stable mode of the critical solution in each case ||l^,|ll|,|| . 
Of the examples considered to date, the evolution of the 
Yang-Mills field is exceptional. Choptuik et al |lj] have 
studied this model, finding two distinct types of phase 
transition depending on the initial field configurations 
they considered. In what they refer to as Type I tran- 
sitions, the black hole formation turns on at finite mass 
and the critical solution is the Bartnik-McKinnon solu- 
tion Q — a regular, static, but unstable, solution to the 
spherically symmetric Einstein- Yang-Mills equations. In 
Type II transitions black-hole formation turns on at in- 
finitesimal mass and the critical behavior is qualitatively 
similar to that found by Choptuik for massless scalar 
fields, except that the scaling exponent is 7 ~ 0.20 and 
the echoing period is A ~ 0.74. An independent con- 
firmation of these results has been provided by Gund- 
lach Q. 

Here we report on a detailed study of critical phe- 
nomenon in the collapse of spherically symmetric config- 
urations of a massive scalar field. The introduction of a 
mass /i destroys the scale invariance of the Einstein-scalar 
field equations. Moreover, the massive scalar-field equa- 
tions admit soliton-like solutions as discussed by Seidel 
and Suen |^ . These observations suggest that the qual- 
itative picture of critical point behavior could differ from 
the massless limit, and might be similar to that found by 
Choptuik et al. jl^ in their study of Einstein- Yang-Mills 
collapse^. Our results show that both Type I and Type II 



*In this case, the SU(2) Yang-Mills charge breaks the scale 
invariance of the field equations. 
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phase transitions occur in the collapse of massive scalar 
fields. Furthermore, we advance a simple criterion to de- 
termine which type of phase transition will be observed 
for a given initial data set. If the radial extent A of the 
initial pulse is greater than the Compton wavelength 
of the scalar field then Type I phase transitions will be ob- 
served. Type II transitions develop from initial data with 
A/i ^ 1. This criterion provides an intuitive, physical ex- 
planation of the observed phenomenology, and clarifies 
the role played by intrinsic scales in critical collapse. 

We write the general, spherically symmetric line- 
element in terms of a retarded time u and a radial co- 
ordinate r, which measures proper area of the 2-spheres, 
as 

ds2 ^ -ggdv^ - 2gdudr + r'^dVi^ , (2) 

where g and 'g are functions of both r and u, and dVl^ 
is the line-element on the unit 2-sphere. We choose to 
normalize u to be proper time at the origin, thus fixing 
(7(0, m) = 1. Imposing regularity of the spacetime at the 
origin requires g(0,u) = 1. The evolution of the massive 
scalar field is governed by the wave equation Ucf) — _ 
0. It is convenient to introduce an auxiliary field h p6| , 
related to </> by 

4> = h = - i dr'h{r',u) , (3) 
^ Jo 

in terms of which the wave equation can easily be written 
in a first order form. Thus, the coupled Einstein-scalar 



field equations are 

{\ng),r = r-\h-h)\ (4) 

{rg),r = .g(l - fi'r't) , (5) 

irh),r ^h, (6) 

h^u - = —ih-h)[g{l - y?r^h^)~g] 

-rgi?h/2 . (7) 



The dynamics is completely encompassed in the last 
equation, which is the wave equation written in terms 
of the new fields h and h. 

The characteristic initial value problem requires only 
the field (p to be supplied on some initial outgoing null 
cone, which we will take, without loss of generality, to 
be at u = 0. We have considered the evolution of three 
different initial data sets as shown in Table |. 

The numerical algorithm used to integrate these equa- 
tions is documented elsewhere. We have followed 
the scheme as outlined by Garfinkle The accuracy 
of the code has been tested previously, where it was used 
to study radiative tails of a massless scalar field propa- 
gating in asymptotically de Sitter spacetimes p^ , and 
was found to be locally second order accurate. 



TABLE I. The three initial data sets considered in our 
evolutions. Only a single parameter is varied when looking 
for a critical point. The types of phase transitions which may 
occur are indicated under Type. 



(j){u — 0, r) Parameters Type 

(i) (j)or^ exp[-{r - rof /a^] a, (j)o I, II 

(ii) 00 {1 - tanh[(r- - ro)/cr]} a, (po I, II 

(iii) <?ior(r + r-o)-"/(l + e'-) a, 9^0 I, II 



It is reasonable to expect that the picture of critical be- 
havior offered by Choptuik is robust against the intro- 
duction of a small scalar-field mass [^-|| . More precisely, 
the evolution of an initial distribution of scalar field will 
differ from the massless evolution only if the character- 
istic length-scale set by the scalar field mass, is 
smaller than the radial extent A of the region in which 
the field is non-zeroQ. This expectation is supported by 
our numerical integrations. We generally observe Type 
II phase transitions when A/i <^ I. Furthermore, we find 
that the scaling relation for the black hole mass is in 
agreement with the massless limit, having 7 « 0.378. 

The new feature, arising due to the presence of the 
mass /i, is the existence of Type I phase transitions — 
phase transitions in which black hole formation turns on 
at finite mass. The critical solutions are soliton stars 
on the unstable branch of the mass versus radius curve 
discussed by Seidel and Suen The mass gap at 

the threshold of black hole formation lies in the range 
0.35 < iiMbh ^ 0.59, the upper limit being set by the 
maximum mass that a soliton star can have. 

For the three initial data sets we examined (see Ta- 
ble I), we have found both Type I and Type II behavior, 
along with evidence that both critical solutions play a 
role when A/i « I. In contrast to the Einstein- Yang- 
Mills system our results suggest that the shape of 
the initial data does not determine the critical point be- 
havior. 

Physically the existence of the different regimes can be 
understood if we recall the two known limits for scalar 
field collapse. In the massless regime A/i ^ I an outward 
pressure is required for the field to bounce back to infin- 
ity, whereas in the adiabatic regime A/i ^ I the collapse 
is pressureless |2^. Hence, by continuity in the space 
of solutions, it seems likely that there can exist config- 
urations, characterized by A/i = C ^ C'll); such that 
the field neither disperses to infinity nor collapses to a 



*For generic initial data it is not possible to define the radial 
extent, however, for the initial data sets (i) and (ii) in Table ^ 
we set A = 2a. 
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singularity. 

Figure |l| shows the field (p at the origin r — during 
the critical phase of a Type I evolution for the Gaus- 
sian initial data in Table | with tq = 5.0, a = 1.2 and 
00 — 6.4746. The solution behaves like a soliton star 
with an effective mass ~ 0.52/i~^ confined within a ra- 
dius 4/i~^ for an amount of time u ~ 100. The an- 
gular frequency of its fundamental mode is ujq w 1.8//, 
the next mode is also apparent at Swq- Superimposed on 
these oscillations is an amplitude modulation with period 

15.7/x~^ which is reflected in the sideband structure of 
the Fourier amplitude. 

Sub-critical evolutions generally settle down to a scalar 
field configuration dominated by a single oscillatory mode 
with angular frequency uj w 1.05 as shown in Fig ^ for 
fi ~ 1.0. Further exploration indicates that the funda- 
mental oscillations have a period given approximately 
by 27r/i^^. Unfortunately, the characteristic evolution 
scheme makes it difficult to follow the evolutions beyond 
u ~ 400. Nevertheless, the trend suggests that the am- 
plitude of the scalar field is decreasing slowly — indeed 
this is precisely the regime in which the stationary phase 
approximation is valid pO| ] . 

When the critical solution corresponds to the 
marginally stable soliton star with effective mass ~ 
0.6/i~^, we have found evidence of further phenomenol- 
ogy. In particular, the solutions may closely approach 
the solitonic conffguration, begin to disperse, but recol- 
lapse to form black holes. This behavior merits further 
investigation, however our numerical scheme is not well 
suited for this purpose. 

As evidence of the observed mass gap, we present in 
Fig. ^ the spectrum of black hole masses near to criti- 
cality for the initial data set (iii) of Table |. The black- 
hole mass at threshold is Mbh ~ O.blfi . This mass 
spectrum is most interesting for what it does not do, 
rather than what it does. We determine that a solution 
contains a black hole if either of the metric functions 
g 01 g exceeds some pre-specified tolerance Gmax any- 
where on a slice of constant u. Suppose this occurs at 
u = ubh[4'o]: where (f>o is the parameter being varied 
in the initial data. The mass of the black hole is then 
Mbh — \''^BH where tbh is the location of the global 
minimum of f{r\(j)(^) = g{r,UBH[(t}o])/ gir,UBH[(f>o]) on 
the slice u — ubh- Notice that Mbh need not depend 
continuously on (po if f{r;(f)o) has more than one local 
minimum. Indeed, as 0o is varied in our simulations 
we sometimes observe discontinuities in the black-hole 
mass when the tolerance Gmax is not large enough; this 
is shown in the inset of Fig. ^. A careful inspection, 
varying both the numerical resolution and the tolerance, 
suggests that the discontinuities are not a real effect in 
the black-hole mass spectrum. (Discontinuities in the 
mass spectrum are also alluded to in p^ ; it would be 
interesting to check if they arise for similar reasons.) In 
contrast, the oscillation imposed on the mass spectrum 



in Fig. ^ is not an artifact of the numerics but is similar 
to the fine structure found by Hod and Piran [ pl[ in the 
Choptuik results. 

To better understand the selection effect between Type 
I and Type II phase transitions we have constructed fam- 
ilies of interpolating solutions [4>o] for several values of 
A. Generally, we find Type I transitions occur when the 
Bondi mass Afeondi of the initial field profile is greater 
than ^ 0.4/i^^ and its radial extent is larger than the 
Compton wavelength of the field, i.e. X/i ^ 1. Figure ^ 
shows the Bondi mass of the initial data at the critical 
point 00 = 00, and the resulting black- hole mass, for the 
initial data sets (i) and (ii) in Table |. 

In conclusion, we find that the presence of a length 
scale changes the nature of critical phenomena in gravi- 
tational collapse of a scalar field. It introduces new phe- 
nomenology which is similar to that discussed by Chop- 
tuik et al jljl . Moreover, it is tempting to speculate that 
unstable, confined solutions will act as critical solutions 
in other matter models. We therefore expect that both 
Type I and Type II phase transitions should occur in the 
gravitational collapse of perfect fiuids (with equations of 
state which allow stationary configurations), and in the 
collapse of charged massive scalar fields [|2| . 
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FIG. 1. (a) The time evolution of as a function of re- 
tarded time u at the origin for a near critical evolution 
with Gaussian initial data (See Table ^ . It clearly exhibits 
an underlying periodic solution with a superimposed ampli- 
tude modulation, (b) The squared amplitude of the discrete 
Fourier transform of (/)(0, it). The fundamental oscillation has 
an angular frequency ujq ~ 1.8/^, and, in agreement with Sei- 
del and Suen |l5| , the next important feature is at three times 
this frequency ui ~ 5.65/i. The sidebands determine the pe- 
riod of the amplitude modulation to be ~ 15.7 fj,~^ . 
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3. The black-hole mass Mbh as a function of 
- (pol for supercritical evolutions with jj. — 1.0. The re- 
sults displayed are for the initial data set (iii) in Table | with 
ro = 2.0 and cr — 10.0. The critical point was determined 
to be (pQ = 3.87245233459 with an initial radial discretization 
Ar = 0.05. The black hole tolerance parameter (see text) was 
set at lO^''. The inset shows results obtained from the same 
evolution, but with a lower black hole tolerance level of 10*. 
For low values of the tolerance the mass spectrum exhibits 
spurious discontinuities. 
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FIG. 2. (a) The time evolution of as a function of re- 
tarded time u at the origin for the same evolution as in Fig |l| 
but after the nearly critical phase which ends at u ~ 100. The 
solution is essentially periodic, with a low frequency ampli- 
tude modulation, (b) The squared amplitude of the discrete 
Fourier transform of the field for it > 100. The fundamental 
oscillation has an angular frequency ujo ~ 1.05/i, which sug- 
gests that it is determined by the scalar field mass (which is 
set to unity). The inset demonstrates that other harmonics 
are not strongly excited in this solution. For other initial data 
sets, we find that the amplitude of the oscillations tends to 
decay slowly, and the solutions are dispersing. 
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FIG. 4. The Bondi mass of the initial scalar field profile, 
in practice m(r, 0) = r(l — g/g)/2 evaluated at the outer 
edge of the grid, and the measured black-hole mass at the 
critical point versus the radial extent A of the initial profile 
for the data sets (i) and (ii) in Table |. Type I transitions are 
evident for Xfi S> 1, and Type II transitions when Xfi <g 1. 
The interface between Type I and Type II behavior is clearly 
visible when -MBondi — 0.4^~^, and Xfj, ~ 1. 
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